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The semiclassical dynamics of a charged particle moving in a two-component plasma is consid- 
ered using a corrected Kelbg pseudopotential. We employ the classical Nevanlinna-type theory of 
frequency moments to determine the velocity and force autocorrelation functions. The constructed 
expressions preserve the exact short and long-time behavior of the autocorrelators. The short-time 
behavior is characterized by two parameters which are expressable through the plasma static correla- 
tion functions. The long-time behavior is determined by the self-diffusion coefficient. The theoretical 
predictions are compared with the results of semiclassical molecular dynamics simulation. 
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I. INTRODUCTION 

The purpose of this paper is the investigation of the dynamics of force on a charged particle in a two component 
plasma. Boercker et al. have shown the effect of ion motion on the spectral line broadening by the surrounding 
plasma . In recent papers it was argued that the microfield dynamics influence the fusion rates || and rates for 
three-body electron-ion recombination Q in dense plasmas. Generally speaking, to calculate the plasma effect on 
rates and spectral line broadening one needs a theory of average forces and microfields, including the resolution in 
space and time. Basic results in this field were obtained by Silin and Rukhadze ||, Klimontovitch ||, Alastuey et. 
al. 0, and Berkovsky et. al. ||. 

The determination of the static distribution of the ion or electron component of the electric microfield is a well 
studied problem (for a review see The corresponding investigations are performed on the basis of the one- 

component plasma (OCP) model. A straightforward generalization of the OCP model is the model of a two-component 
plasma (TCP), consisting of electrons and ions. In a recent paper JlCfl the probability distribution for the electric 
microfield at a charged point has been studied. It was shown that the two-component plasma microfield distribution 
shows a larger probability of high microfield values than the corresponding distribution of the OCP model. 

The dynamics of the electric microfield is a less understood problem than that of the static microfield distribution 
even for the case of an OCP. Recently some progress has been made for both the case of electric field dynamics at a 
neutral point |f^, p] , ^2[ and the dynamics of force on a charged impurity ion in an OCP Q . 

This paper is aimed to extend the studies of electric microfield dynamics in OCP to the case of an equilibrium two 
component plasma. For simplicity we consider a two-component plasma which is anti-symmetrical with respect to 
the charges (e_ = — e + ) and therefore symmetrical with respect to the densities (n^ = n e ). To simplify the numeric 
investigations we simulated a mass symmetric (nonrelativistic) electron-positron plasma with m = m t = m e . The 
theoretical investigations are carried out for arbitrary electron-ion mass ratios. 

1 Dedicated to the 75th birthday of Youri L. Klimontovich 
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In this paper we will study the dependence of the force dynamics on the coupling constant T = e 2 /ksTa of the 
plasma, where T is the temperature, and a — (3/47rn e ) 1 / 3 is the average distance between the electrons. Coupled 
plasmas with a plasma parameter of the order or greater than unity are important objects in nature, labora tory 
experiments, and in technology |]l3|-|l6|. Recent lasers allow to create a coupled plasma within femtoseconds [ ^7[ . 
Laser generated plasmas are nonequlibrium plasmas with an initial electron temperature much higher than the ion 
temperature. However, in this paper we restrict our considerations to the model object of an equilibrium two- 
component plasma (TCP). 

Several investigations were devoted to the simulation of equilibrium two-component plasmas Being interested in 
quasi-classical methods we mention explicitely the quasi-classical simulations of two-component plasmas performed 
by Norman and by Hansen f3~8|J]~s| ] . 

In this paper the free charges (electron and ions) are simulated by a semi-classical dynamics based on effective 
potentials. The idea of the semi-classical method is to incorporate quantum- mechanical effects (in particular the 
Heisenberg and the Pauli principle) by appropriate potentials. This method was pioneered by Kelbg, Deutsch and 
others Certainly, such a quasi-classical approach has several limits. For the calculation of a standard macro- 

scopic property as the microfield dynamics which has a well defined classical limit the semi-classical approach may be 
very useful. The advantage of such an approach is the relative simplicity of the algorithm. 



II. THE SLATER SUM AND THE SEMICLASSICAL MODEL 



A familiar derivation of effective potentials describing quantum effects is based on the Slater sums which are defined 
by the N - particle wave functions, 

S(ri, . . . ,r N ) = constj^exp (-/? E n ) |*„ (n, . . . ,r N )\ 2 , (1) 

where E n and ^ n are the energy levels and corresponding wave functions of the ensemble of N particles with coor- 
dinates r\, . . . , rjsf. Here we consider a two-component plasma consisting of N e electrons with mass m e and Ni = N e 
ions with mass m,. The properties of the Slater sums for Coulombic systems were studied in detail by several authors 
|^3|,p2[ . Choosing the effective potential 

U^ N \r 1 ,...,r N ) = -k B T\nS(r 1 ,...,r N ) . (2) 



we may calculate the correct thermodynamic functions of the original quantum system |13j , |22| , |18| ] from the thermody- 
namic functions of a classical reference system. 

The Slater sum may be considered as an analogue of the classical Boltzmann factor. Therefore it is straightforward to 
use the Slater sum for the definition of an effective potential. The only modification in comparison with classical theory 
is the appearance of many-particle interactions. If the system is not to dense (i.e., n e A 3 e < 1, A e = h/y/2m e k B T) one 
may neglect the contributions of higher order many-particle interactions. In this case one writes approximately, 

uW(r 1 ,...,r N )*Y l u i j(r i ,r j ) , (3) 

i<j 

where the effective two-particle potential u ab is defined by the two-particle Slater sum, 

S ( ab (T) = exp (-/?u tt6 (r)) - const.^' exp {-0E a ) | * Q | 2 . (4) 

a 

Here f a and E a denote the wave functions and energy levels of the pair ab, respectively. The prime at the 
summation sign indicates that the contribution of the bound states (which is not be considered here) has to be 
omitted. 

Principal it is possible to calculate the Slater sum for a pair of particles directly from the known two-particle 
Coulomb wavefunctions. To simplify the simulations it is better to have an analytic expression for the potential. A 
possible candidate is the so called Kelbg potential obtained by a perturbational expansion It reads [|o| 

u ab (r) = —F(r/X ab ) , (5) 
r 

where X a b — h/y/2m a bkBT is De Broglie wave length of relative motion, = m^ 1 + m^ 1 , a = e,i. In Eq.([|) 

F{x) = 1 - exp (-x 2 ) + Vkx (1 - erf (x)) . (6) 
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Another analytic approximation for the exact two-particle effective potential is the expression derived by Deutsch 
which was used in the simulations by Hansen and McDonald |ll| . 

The Kelbg potential is a good approximation for the two-particle Slater sum in the case of small parameters 
tlab = — (e a e(,)/(fcsTA a h) if the interparticle distance r is sufficiently large. However, at small interparticle distances 
it shows a deviation from the exact value of — k B T ■ \n(S ab (r = 0)). In order to describe the right behavior also at 
small distances it is better to use a corrected Kelbg potential defined by 0] 



u a b(r) = {e a e b /r) ■ \F(r/\ ab ) - r^^A ab {^ ab ) exp (-(r/A afc ) 2 ) 1 



(7) 



In Eq. (^) the coefficient A ab {T) is adapted in such a way that S ab {r = 0) and his first derivative S' ab (r — 0) have 
the exact value corresponding to the two-particle wave functions of the free states JlJ,|24|,^3|] . The corresponding 
coefficients for the clctron-clcctron and for the clcctron-ion interaction read 



A ee = V^l^ee 
A e i = — 



hi 



hi 



dy y exp (-y 2 ) 
exp(7r|£ ee |/y) - I 



yfttl ( C(3) + i C(5)eL ) + v< e 



dyyexp (-y 2 ) 
1 - exp (-ttCm/v) 



We mention that in the region of high temperatures 



T r = T/Ti = (2k B Th 2 /m le e A ) > 0.3 



(8) 
(9) 

(10) 



the Kelbg potential (A ab = 0) almost coincide with the corrected Kelbg potential Eq. (|7|). In the region of intermediate 
temperatures 0.1 < T r < 0.3 the Kelbg potential does not give a correct description of the two-particle Slater sum 
at short distances. Instead we may use the corrected Kelbg-potential Eq.(Q) to get an appropriate approximation for 
the Slater sum at arbitrary distances. 

The effective potentials derived from perturbation theory do not include bound state effects. The other limiting 
case of large g a fc or small temperature T r < 0.1, where bound states are of importance, can be treated by another 
approach p2fl. Here a transition to the chemical picture is made, i.e. bound and free states have to be separated. 

In the present work we are interested in the regime of intermediate temperatures. In this regime the simulations of 
the dynamics may be performed with the potential Eq.(^). 



III. FORCE-FORCE AUTOCORRELATION FUNCTION 



The system under consideration is a two-component plasma consisting of electrons and ions which is described by 
the semiclassical model introduced in Sec [n]. Let us choose the position of one of the charged particles (for example 
an electron) as a reference point. Hereafter we call this particle the first one. The semiclassical force acting on the 
first particle equals 



N 

-Ax^uij-Cn - rj ) (li) 

J=2 



ij being the effective pair potential between the iih and jth particles, defined in Eq. (0). 
Define now two functions characterizing the dynamics of the first particle. The first one 



at) = < jKg) > (12) 

< V 2 > 

is the velocity- velocity autocorrelation function (velocity acf), the second function 

C(t) = <m_m> (13) 

is the force- force autocorrelation function (force acf). In the above equations the brackets < . . . > denote averaging 
over the equilibrium ensemble of the semiclassical system. The velocity acf is formally a function expressing the single 
particle properties. However, it is connected with the force acf which involves the collective properties by the relation 
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d 2 c{t) 
at 2 



+ <4D(t)=0, (14) 



where u\ — < F 2 >/3mk B T. 

Define the one-side Fourier transform of the velocity and force acf, 



C(oj) = I dte lu)t C{t) , D(uj) = / dte lu)t D{t) . (15) 
Jo Jo 

The Fourier transform of Eq. (|l4|) reads 

» , , lo 2 C(lo) — iuj . . 

D(u) = 4 . (16) 

In order to construct the both autocorrelation functions it is useful to consider the frequency moments of the real 
part of the velocity acf Fourier transform 

i r 00 

M « = 7T / u n C r (uj)e- lult duj , n = 0, 1, 2, . . . . (17) 



2tt . 

The zeroth moment is the initial value of the velocity acf, 

Mo = C(0) = 1 . (18) 

Due to the parity of the function CV(lj), all moments with odd numbers are equal to zero. 
The second moment is expressable through the initial value of the force acf, 

i r°° 

M 2 = — / uj 2 C r (uj)e- luJt duj = uj 2 D(0) = uj\ . (19) 
The fourth moment includes the correlation function of the time derivative of the force, 

2tt 



M 4 = -!- / ^Cri^e-^dtu = uluo 2 2 , (20) 



where we have introduced the magnitude oj 2 — < F 2 >/< F 2 >. 

The Nevanlinna formula of the classical theory of moments [ ^5|j26| expresses the velocity acf Fourier transform 



1 f°° CV(w) -A/ \ E n+1 (z) + q n (z)E n (z) 

-duj = -iC{z) = — — -j— (21) 



7T z-lj ' D n+ i(z) + q n (z)D n (z) 

in terms of a function q n = q n (z) analytic in the upper half-plane Imz > and having a positive imaginary part there 
Im<| n (w + irf) > 0, rj > 0, it also should satisfy the limiting condition: (q n (z)/z) — > as z — * oo within the sector 
9 < arg(z) < 7r — 6. In Eq.(^l|) we have employed the Kramers-Kronig relation connecting the real and imaginary 
part of C(lu). The polynomials D n (and E n ) can be found in terms of the first 2n moments as a result of the Schmidt 
orthogonalization procedure. The first orthogonal polynomials read 

Di=z, D 2 = z 2 -uj 2 , D 3 = z(z 2 - uj 2 ) , (22) 
E x =l, E 2 = z, E 3 = z 2 + uj 2 - uj 2 ). (23) 

Consider first the approximation n = 1 leading to the correct frequency moments Mq and M 2 - Using the Nevanlinna 
formula and Eq. ( |l6| ) we obtain 

The physical meaning of the function q\(z) is that of a memory function j^] since the inverse Fourier transform of Eq. 

*<j® +u %C{t)+ r t dsqi (t-s)^-=Q. (25) 
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We have no phenomenological basis for the choice of that function qi(z) which would provide the exact expression for 
C(z) and D(z). A simple approximation is to put the function qi(z) equal to its static value 

qi (z) = qi (0) = iv (26) 



and Eq. (25) simplifies to the equation of a damped oscillator with frequency uj\ and damping constant v. 

The static value q\{z = 0) is connected with the self-diffusion coefficient D. The latter is defined by the time 
integral of the velocity acf 

D =ir* c7(t)= ^ (o) ' (28) 



where (3 = l/(fcsT) and m\ is the mass of the first particle. With the use of Eqs. ( |28|) and (26) we obtain from Eq. 
([H]) that v = LujPmxD. 

The inverse Fourier transform of Eq. (|l|) with the static approximation Eq. (|2^) expresses the velocity and force 
acf's as a linear combination of two exponential functions exp(zii) and exp^i), where zi/ 2 = ± \/ v 2 — 4w 2 /2. 

Within this approximation we may distinguish between two regimes. In the first regime - the "diffusion-regime" - one 
deals with a large diffusion constant. As a result v = (3ra\Dio\ > 2ui\ and Eq. ( pTj ) is the equation of an overdamped 
oscillator. In this regime the velocity autocorrelation function goes monotoneously to zero. With decreasing diffusion 
constant the damping constant v becomes smaller. At certain thermodynamical conditions just the opposite inequality 
v < 2uj\ holds. This corresponds to an "oscillatory-regime" and at least one of the autocorrelation functions should 
show an oscillatory behavior. The existence of the two regimes have been established for the case of an OCP J8| and 
has been confirmed by our molecular-dynamics simulation for the case of a TCP. To obtain not only a qualitative but 
also a quantitative correspondence with the results of MD simulations one has to go beyond the simple approximation 
n = 1 in the Nevanlinna formula Eq. (pl|). 

Consider therefore the case n = 2 in Eq. (^lj). Then the autocorrelation functions are expressed via the function 
q 2 (z) as 

C(2) = %(z 2 -i 2 ) + g 2 (^-^) ' D{z) = t z(z 2 -u;l)+ q2 (z 2 ~Lul)- (M) 

Eq. ( |29| ) reproduces the exact freqency moments from M$ up to M4. For the function q 2 (z) we choose again a static 
approximation 

q 2 (z) = q 2 {0) = ih, (30) 

where h has to be taken from the relation 

h = - l\ /(3m x D (31) 



in order to obtain the exact low frequency value C(0) given by Eq. 

From Eq. (29) we find that the autocorrelation functions are now given by the linear combination of three expo- 
nentials, 



3 3 



C(t) = ) Qe'"'* , D(t) = > d.^ . (32) 



= 1 i=l 



The complex frequencies il^ are the poles of the expressions Eq. ([29|) . They are defined as the solutions of the cubic 
equation, 

Vl{n 2 - ujf) + ih{tt 2 - u\) =0. (33) 
The coefficients Ci (di) characterizes the strength of the zth mode, 



5 



C t = j£di, i = 1,2,3, (34) 

d 1 =i{h + in 1 )Q 1 (Q 2 -^3)/N (35) 

d 2 = i(h + in 2 )n 2 {n 3 -n 1 )/N (36) 

d 3 = i(h + ^3)^3(^1 (37) 

J\r=(ni-n 2 )(n3-ni)(fi2-fi3). (38) 



Equations ( p2[ ) constitute the basic approximation of our paper. The frequencies f2j and the coefficients C{ (or 
di, respectively) are expressed by three parameters - the diffusion constant D, and the frequencies lo\ and uj%. The 
constructed autocorrelation functions satisfy the following conditions: (i) the exact short time behavior for the velocity 
acf is reproduced to the orders t 2 and t 4 , (ii) the short time behavior of the force acf is reproduced to the order t 2 , 
(iii) the long time behavior of the velocity acf generates the exact diffusion constant, and (iv) the connection between 
the velocity and force acf 's Eq. ( |l4"| ) is satisfied. 

The parameters D, lo\ and oj 2 may be calculated by another approximations. The both frequencies u)\ and u> 2 
are expressable via the partial correlation functions of our semiclassical system. The parameter uj\ is given by the 
electron-ion and electron-electron partial pair correlation functions. To calculate the frequency lo 2 one needs the 
knowledge of the partial ternary distribution functions. The diffusion constant may be obtained from kinetic theory. 
In contrast to the case of an OCP Q the parameters to be calculated are very sensitive to the approximations used 
to calculate the static distribution functions. Therefore in this paper we take the "input" parameters directly from 
the computer simulations. 

To check the quality of the predictions from our approximation we have performed molecular dynamics simulations 
for comparison. The equations of motions obtained with the effective potential Eq. (^) were integrated numerically for 
the case of equal masses m e = rrii using the leap-frog variant of Ver let's algorithm. The simulations were performed 
for 128 electrons and 128 positrons moving in a cubic box with periodic boundary conditions. In the investigated 
range of plasma parameters (T = 30 000 K, the coupling parameter has varied from T = 0.2 up to T = 3) the size 
of the simulation box was significantly greater than the Debye radius. Therefore the long-range Coulomb forces are 
screened inside each box and there was no need to use the Ewald summation instead the simple periodic boundary 
conditions. The thermal equilibrium in the system was established (and maintained) by a Langevin source. Such 
simulations has been recently used to obtain the static distribution of the electric microfield at a charged particle [[l0| . 
In this paper we extract the velocity and force autocorrelation functions as the main characteristics of the microfield 
dynamics. 

TABLE I. The T dependence of the parameters lo 2 and D. uj\ and io 2 are given in units of electron plasma 
frequency aj pe = y/ 4irn e e 2 /m e , D is given in units of 1 / (m e aj pe (3) 



r 


wi 


LU2 


D 


0.2 


0.84 


13.6 


10.3 


1.5 


0.88 


3.3 


4.41 


3.0 


0.61 


2.1 


5.75 



In Figs. 1-3 we present the results of the MD data. The simulation results are compared with our analytical 
approximation Eqs. ((2^)- The three input parameters for the analytical approximation are taken from the MD 
simulations. The diffusion constant is obtained from the time integral of the velocity acf (Eq. (28)). Since the 
velocity acf is a slowly decaying function it requires a long simulation time to extract the diffusion constant. For our 
model system with equal electron and ion masses it is possible to perform the necessary simulations. The frequency 
u> 2 has been taken from the exact short time behavior of the force acf D(t) = 1 — uj\t 2 /2. Finally the frequency u>\ was 
choosen to fit the model to the data. In Table I we show the parameters lo\ , oj 2 and D for three coupling parameters 
r considered in this paper. 

Except the case of the force acf at T = 0.2 there is a good overall agreement between the theoretical approximations 
and the MD data. We believe that the strong deviation of the MD data from the theoretical predictions for T — 0.2 is 
a numerical artefact due to the poor statistics in the weak coupling case. From the figures we see that with increasing 
plasma parameter T the dynamics of the charged particles switches from the diffusion-like regime at T = 0.2 to the 
oscillator-like motion at T = 3.0. The value T = 1.5 may be considered as a critical value separating the both regimes. 
We may also see from the figures that the oscillator-like motion is more pronounced for the force acf. 

At still higher densities (r > 3 at T = 30 000 K) the semi-classical approach employed in this paper fails to describe 
the quantum two-component plasma properly. 
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IV. CONCLUSIONS 



The electric microfield dynamics at a charged particle in a two-component plasma has been studied. The quantum 
plasma has been modeled by a semiclassical system with effective potentials. The effective potential was choosen to 
describe the nondegenrate limit of the quantum system appropriately. We have investigated the velocity and force 
aef's of the semiclassical system. The starting point for the theoretical analysis was the exact expression of the 
autocorrelation functions through the Nevanlinna formula Eq. (|29|), satisfying three sum rules for the velocity acf. 
The approximation Eq. ( p0| ) together with Eq. ( |3l| ) expresses the velocity acf in terms of three parameters. Two of 
them - oji and lo 2 - describe the exact short time behavior of the velocity acf up to the order i 4 , the third parameter, 
the self-diffusion constant D is related to the time integral of the velocity acf. Since the force acf can be obtained 
from the velocity acf by a second time derivative the force acf is expressed through the same three parameters. The 
general picture is as follows. At weak coupling the diffusion of the charged particle dominates the collective plasma 
oscillations and the particle motion is diffusion-like. In this regime the velocity acf decays exponentially with a decay 
rate 1/D (time in units of the inverse electron plasma frequency u> pe ). The force acf has a positive decay at short 
times (decay rate to\D) and a negative decay at long times (with the rate At strong coupling the diffusion is 

supressed and a weakly damped oscillatory behavior for the force acf developes. The theoretical predictions has been 
compared with molecular dynamics simulations data. There is an overall agreement of the force dynamics obtained 
by the analytical approximation with the MD data. 

Finally, wc mention that there is no one to one correspondence of the semiclassical autocorrelation functions with the 
corresponding characteristics of the quantum system. Nevertheless, we suspect that the semiclassical force dynamics 
considered in this paper at least qualitatively reproduces the electric microfield dynamics of the quantum system. 
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FIGURE CAPTIONS 



(Figure 1) Time dependence of velocity acf C(t) and force acf D(t) at T = 0.2. Time is in units of inverse electron 
plasmafrequency c^ 1 . Solid lines: present theoretical approximation; Points: results of molecular-dynamics 
simulations. 

(Figure 2) Same as in Fig. 1 at T = 1.5. 

(Figure 3) Same as in Fig. 1 at T = 3.0. 
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